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Summary. We consider spacetime to be a 4-dimensional differentiable manifold 
that can be split locally into time and space. No metric, no linear connection are 
assumed. Matter is described by classical fields/fluids. We distinguish electrically 
charged from neutral matter. Electric charge and magnetic flux are postulated to 
be conserved. As a consequence, the inhomogeneous and the homogeneous Maxwell 
equations emerge expressed in terms of the excitation H — (Ti, T>) and the field 
strength F — (E, B), respectively. H and F are assumed to fulfill a local and linear 
"spacetime relation" with 36 constitutive functions. The propagation of electromag- 
netic waves is considered under such circumstances in the geometric optics limit. 
We forbid birefringence in vacuum and find the light cone including its Lorentzian 
signature. Thus the conformally invariant part of the metric is recovered. If one sets 
a scale, one finds the pseudo-Riemannian metric of spacetime. 

1 Introduction 

The neutrinos, in the standard model of elementary particle physics, are as- 
sumed to be massless. By the discovery of the neutrino oscillations, this as- 
sumption became invalidated. The neutrinos are massive, even though they 
carry, as compared to the electron, only very small masses. Then the photon 
is left as the only known massless and free elementary particle. The gluons 
do not qualify in this context since they are confined and cannot exist as free 
particles under normal circumstances. 

Consequently, the photon is the only particle that is directly related to the 
light cone gijdx l ®dx : > — and that can be used for an operational definition of 
the light cone; here is the metric of spacetime, dx % a coordinate differential, 
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and i, j = 0,1, 2, 3. We are back — as the name light cone suggests anyway 
— to an electromagnetic view of the light cone. Speaking in the framework 
of classical optics, the light ray would then be the elementary object with the 
help of which one can span the light cone. We take "light ray" as synonymous 
for radar signals, laser beams, or electromagnetic rays of other wavelengths. 
It is understood that classical optics is a limiting case, for short wavelengths, 
of classical Maxwell-Lorentz electrodynamics. 

In other words, if we assume the framework of Maxwell-Lorentz electrody- 
namics, we can derive, in the geometric optics limit, light rays and thus the 
light cone, see Perlick [48] and the literature given there. However, the formal- 
ism of Maxwell-Lorentz electrodynamics is interwoven with the Riemannian 
metric g^ in a nontrivial way. Accordingly, in the way sketched, one can never 
hope to find a real derivation of the light cone. 

Therefore, we start from the premetric form of electrodynamics, that is, a 
metric of spacetime is not assumed. Nevertheless, we can derive the generally 
covariant Maxwell equations, expressed in terms of the excitation H = (H,T>) 
and the field strength F = (E,B), from the conservation laws of electric 
charge and magnetic flux. We assume a local and linear spacetime relation 
between H and F. Then we can solve the Maxwell equations. In particular, 
we can study the propagation of electromagnetic waves, and we can consider 
the geometrical optics limit. In this way, we derive the light rays that are 
spanning the light cone. In general, we find a quartic wave covector surface 
(similar as in a crystal) that only reduces to the pseudo-Ricmannian light 
cone of general relativity if we forbid birefringence (double refraction) in vac- 
uum. Hence, in the framework of premetric electrodynamics, the local and 
linear spacetime relation, together with a ban on birefringence in vacuum, 
yields the pseudo-Ricmannian light cone of general relativity. Accordingly, 
the geometrical structure of a Riemannian spacetime is derived from purely 
electromagnetic data. We consider that as our contribution to the Einstein 
year 2005, and we hope that going beyond the geometrical optics limit will 
yield better insight into the geometry of spacetime. 

The axiomatic scheme that we are going to present here is already con- 
tained in our book [20] where also references to earlier work and more details 
can be found. In the meantime we learnt from the literature that appeared 
since 2003 (see, e.g., Delphenich [8, 9], Itin [24], Kaiser [26], Kichn [29], and 
Lindell & Sihvola [35, 37]) and improved our derivation of the light cone, 
simplified it, made it more transparent (see, e.g., [21, 23, 33, 45, 46]). The 
formalism and the conventions we take from [20]. 

2 Spacetime 

In our approach, we start from a 4-dimensional spacetime manifold that is 
just a continuum which can be decomposed locally in (1-dimcnsional) time 
and (3-dimensional) space. It carries no metric and it carries no (linear or 
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affinc) connection. As such it is inhomogeneous. It doesn't make sense to 
assume that a vector field is constant in this continuum. Only the constancy 
of a scalar field is uniquely defined. Also a measurement of temporal or spatial 
intervals is still not defined since a metric is not yet available. 

In technical terms, the spacetime is a 4-dimensional connected, Hausdorff, 
paracompact, and oriented differential manifold. On such a manifold, we as- 
sume the existence of a foliation: The spacetime can be decomposed locally 
into three-dimensional folios labeled consecutively by a monotonically increas- 
ing "prototime" parameter a, see Fig. 1. A vector field n, transverse to the 




Fig. 1. Local spacetime foliation, see [20]. 



foliation, is normalized by n\da — C n a — 1. Accordingly, we find for the 
dimensions [n] = [a]^ 1 =t~ 1 , where t denotes the dimension of time. 

We can decompose any exterior form 'P in "time" and "space" pieces. The 
part longitudinal to the vector n reads 

J ~^ r := da A £j_ , &± := , (1) 

the part transversal to the vector n 

£:= (1 - ± )f = n\{daA&), n\ £ = . (2) 

Putting these two parts together, we have the space-time decomposition 

if = + £ = da A + £ , (3) 

with the absolute dimensions [£J = W] t' 1 and [£] = [&]. 

The 3-dimensional exterior derivative is defined by d := n\(da A d). We 
can use the notion of the Lie derivative of a p-form along a vector field £, 
i.e., C^P := £,\d<P + d(£\\t r ), and can introduce the derivative of a transversal 
field £ with respect to prototime as 



£:=£„£. 



(4) 
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3 Matter — electrically charged and neutral 

We assume that spacetime is "populated" with classical matter, either de- 
scribed by fields and/or by fluids. In between the agglomerations of matter, 
there may also exist vacuum. 

Matter is divided into electrically charged and neutral matter. Turning to 
the physics of the former, we assume that on the spatial folios of the manifold 
we can determine an electric charge Q as a 3-dimensional integral over a charge 
density and a magnetic flux # as a 2-dimensional integral over a flux density. 

This is at the bottom of classical electrodynamics: Spacetime is filled with 
matter that is characterized by charge Q and by magnetic flux <P. For neutral 
matter both vanish. The absolute dimension of charge will be denoted by q, 
that of magnetic flux by <f) = [action/charge] = h/q, with h as the dimension 
of action. 



4 Electric charge conservation 

One can catch single electrons and single protons in traps and can count them 
individually. Thus, the electric charge conservation is a fundamental law of 
nature, valid in macro- as well as in micro-physics. 3 Accordingly, it is justified 
to introduce the absolute dimension of charge q as a new and independent 
concept. 

Let us define, in 4-dimcnsional spacetime, the electrical current 3-form </, 
with dimension [J] = q. Its integral over an arbitrary 3-dimcnsional spacetime 
domain yields the total charge contained therein: Q = J n ^ J. Accordingly, the 
local form of charge conservation (Axiom 1) reads: 

dJ = Q. (5) 

This law is metric-independent since it is based on a counting procedure for 
the elementary charges. Using a foliation of spacetime, we can decompose the 
current J into the 2-form of the electric current density j and the 3-form p of 
the electric charge density: 

J = -j Ada + p. (6) 

Then (5) can be rewritten as the continuity equation: 

p + dj = 0. (7) 

Both versions of charge conservation, eqs.(5) and (7), can also be formulated 
in an integral form. 

3 Lammerzahl, Macias, and Miiller [34] proposed an extension of Maxwell's equa- 
tions that violates electric charge conservation. Such a model can be used as a test 
theory for experiments that check the validity of charge conservation, and it allows 
to give a numerical bound. 
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5 Charge active: excitation 



Electric charge was postulated to be conserved in all regions of spacetime. If 
spacetime is topologically sufficiently trivial, we find, as consequence of (5), 
that J has to be exact: 

J = dH. (8) 

This is the inhomogeneous Maxwell equation in its premetric form. The elec- 
tromagnetic excitation 2-form H, with [H] = [J] = q, is measurable with the 
help of ideal conductors and superconductors and thus has a direct operational 
significance. 

By decomposing H into time and space, we obtain the electric excita- 
tion 2-form V (historical name: "dielectric displacement") and the magnetic 
excitation 1-form H ("magnetic field"): 

H = -H A da + V . (9) 

Substituting (9) into (8), we recover the pair of the 3-dimensional inhomoge- 
neous Maxwell equations 

dH = J { . * V=P > (10) 
\-V + dH =j. V ' 



6 Charge passive: field strength 

With the derivation of the inhomogeneous Maxwell equations the information 
contained in Axiom 1 is exhausted. As is evident from the Coulomb-Gauss law 
dT> = p, it is the active character of p that plays a role in this inhomogeneous 
Maxwell equation: The charge density p is the source of V (and, analogously, 
the current density j that of Ti). 

Since we search for new input, it is near at hand to turn to the passive 
character of charge, that is, to wonder what happens when a test charge is 
put in an electromagnetic field. In the purely electric case with a test charge 
e, we have 

F a ^eE a , (11) 

with F a and E a as components of the covectors of force and electric field 
strength, respectively. The simplest relativistic generalization for defining the 
electromagnetic field is then of the type 

force density ~ field strength x charge current density . (12) 

Accordingly, with the force density covector (or 1-form) f a , we can formu- 
late Axiom 2 as 

f a = (e a \F) A J . (13) 
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Fig. 2. Faraday- Schouten pictograms of the electromagnetic field, see [20]. The 
electric excitation V is a twisted 2-form, the magnetic excitation TC a twisted 1- 
form. The electric field strength E is a 1-form and the magnetic field strength B a 
2-form, both without twist. 

Here e a is a local frame, with a = 0, 1, 2, 3. Axiom 2 provides an operational 
definition of the electromagnetic field strength 2-form F, the absolute dimen- 
sion of which turns out to be [F] = h/q. Its 1 + 3 decomposition 



introduces the electric field strength 1-form E and the magnetic field strength 
2-form B, see Fig. 2. If we substitute (14) and (6) into (13), we recover, for 
a = 1, 2, 3, the Lorentz force density. 



7 Magnetic flux conservation 

The field strength F, as a 2-form, can be integrated over a 2-dimensional area 
fli in 4-dimensional spacetime. This yields the total magnetic flux piercing 
through this area: <P — F. In close analogy to electric charge conservation, 
we assume that also the flux is conserved. Then, in local form, magnetic flux 
conservation (Axiom 3) reads 4 



4 One can give up magnetic flux conservation by introducing magnetic monopoles 
according to dF = J mag n . In premetric electrodynamics this has been done by Edelen 
[11], Kaiser [26], and by us [21]. However, then one has to change Axiom 2, too, and 
the Lorentz force density picks up an additional term — (e a J_ff)A Jmagn- This destroys 
Axiom 2 as an operational procedure for defining F. Moreover, magnetic charges 
have never been found. 



F = E Ada + B 



(14) 
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dF = 



dB = 0, 
B + dE = 0. 



(15) 



The Faraday induction law and the sourcelessness of B are the two conse- 
quences of dF = 0. In this sense, Axiom 3 has a firm experimental underpin- 
ning. 

8 Premetric electrodynamics 

...is meant to be the "naked" or "featureless" spacetime manifold, with- 
out metric and without connection, together with the Maxwell equations 
dH = J, dF = 0, the Lorentz force formula, and the electromagnetic energy- 
momentum current to be discussed below, see (20). We stress that the Poincare 
group and special relativity have nothing to do with the foundations of elec- 
trodynamics as understood here in the sense of the decisive importance of the 
underlying generally covariant conservation laws of charge (Axiom 1) and flux 
(Axiom 3). Historically, special relativity emerged in the context of an analysis 
of the electrodynamics of moving bodies [13, 39], but within the last 100 years 
classical electrodynamics had a development of its own and its structure is now 
much better understood than it was 100 years ago. Diffeomorphism invariance 
was recognized to be of overwhelming importance. Poincare invariance turned 
out to play a secondary role only. 

Of course, premetric electrodynamics so far does not represent a complete 
physical theory. The excitation H does not yet communicate with the field 
strength F. Only by specifying a "spacetime" relation between H and F (the 
constitutive law of the spacetime manifold) , only thereby we recover — un- 
der suitable conditions — our normal Riemannian or Minkowskian spacetime 
which we seem to live in. In this sense, a realistic spacetime — and thus an 
appropriate geometry thereof — emerges only by specifying additionally a 
suitable spacetime relation on the featureless spacetime. 

As explained, Axiom 1, Axiom 2, Axiom 3, together with Axiom 4 on 
energy-momentum, constitute premetric electrodynamics. Let us display the 
first three axions here again, but now Axiom 1 and Axiom 3 in in the more 
general integral version. For any submanifolds C3 and C2 that are closed, i.e., 
dCz = and dCi = 0, the axioms read 



By de Rham's theorem we find the corresponding differential versions 






(16) 



dJ = 



f a = (e«J F) A J , 



dF = 



(17) 
(18) 



J = dH, 



F = dA. 
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Fig. 3. Different aspects of the electromagnetic field. The four quantities Ti, T>, E, B 
constitute the electromagnetic field. The excitations Ti, T> are extensive quantities 
(how much?), the field strengths E, B intensive quantities (how strong?). 



The physical interpretation of the quantities involved is revealed via their 
(1 + 3)-decompositions (6), (9), (14), and A = —Lpda + A, see Fig. 3. 

Let us now turn to the energy-momentum question. Using the properties 
of the exterior differential, we can rewrite the Lorentz force density (13) as 

f a = (e a \F)AJ = d k S a + X a . (19) 

Here the kinematic energy-momentum 3-form of the electromagnetic field, a 
central result in the premetric electrodynamics, reads (Axiom 4) 

k £ a :=^[FA (e Q J H ) - H A (e Q J F)} . (20) 

The remaining force density 4-form turns out to be 

X a :=-^(FAC ea H-H AC ea F) . (21) 

The absolute dimension of k S a and of X a is h/£, where i denotes the dimen- 
sion of length. [Provided, additionally, a linear connection is given with the 
covariant differential D, then 

f a = D k S a + X a , (22) 
with the new supplementary force density 

X a = ^ {H A h Ea F — F A L eQ H) , (23) 

which contains the covariant Lie derivative. In general relativity theory, X a 
eventually vanishes for the standard Maxwell-Lorentz electrodynamics.] 
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9 The excitation is local and linear in the field strength 

The system of the Maxwell equations dH = J , dF = is apparently under- 
determined. It gets predictive power only when we supplement it with a space- 
time (or constitutive) relation between the excitation and the field strength. 
As Axiom 5, we postulate a general local and linear spacetime relation 

H = k(F), H ij = ^ Kij kl F kl . (24) 

Here excitation and field strength decompose according to H = Hij dx l f\dx : >/2 
and F = F^ dx % A dx j /2, respectively. The constitutive tensor n, as 4th rank 
tensor with 36 independent components, has to be space and time dependent 
since constant components would not have a generally covariant meaning on 
the naked spacetime manifold we consider. 

Let us decompose Ki kl into irreducible pieces. In the premetric framework 
we can only perform a contraction. A first contraction yields 

n k :— Ku kl (16 independent functions) , (25) 

a second one 

k := Kk k = Kki H (1 pseudo-scalar function) . (26) 
Then, introducing the traceless piece 

^ k := K k - - nS k (15 functions) , (27) 

we can rewrite the original constitutive tensor as 

K .M _ (1) ki , (2) kl , (3) kl 

= ^n ii kl +2^1**2 + ±k5{$]. (28) 
The skewon and the axion fields are conventionally defined by 

^ = -\^, a =^ K - ( 29 ) 

Substituting (28) into (24) and using (29), we obtain the spacetime relation 
explicitly: 

Ha = \ ^ Kij kl F kl + 2 fS { *F j]k + a F iS . (30) 

The principal (or the metric-dilaton) part ^Kij kl of the constitutive tensor 
with 20 independent components will eventually be expressed in terms of the 
metric (thereby cutting the 20 components in half). [In standard Maxwcll- 
Lorentz electrodynamics 
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W Kij kl =X ^e ijmn9 mk g nl , fi k = , a = 0.] (31) 

The principal part W Hij kl must be non- vanishing in order to allow for elec- 
tromagnetic wave propagation in the geometrical optics limit, see the next 
section. The skewon part fi k with its 15 components was proposed by us. 
We put forward the hypothesis that such a field exists in nature. Finally, the 
axion part a had already been introduced in elementary particle physics in 
a different connection but with the same result for electrodynamics, see, e.g., 
Wilczek's axion electrodynamics [65] and the references given there. 

The spacetime relation we are discussing here is the constitutive relation 
for spacetime, i.e., for the vacuum. However, one has analogous structures for 
a medium described by a local and linear constitutive law. The skewon piece 
in this framework corresponds to chiral properties of the medium inducing op- 
tical activity, see Lindell et al. [38] , whereas the concept of an axion piece has 
been introduced by Tellegen [57, 58] for a general medium, by Dzyaloshinskii 
[10] specifically for Cr 2 03, and by Lindell & Sihvola [37] in the form of the 
so-called perfect electromagnetic conductor (PEMC). Recently, Lindell [36] 
discussed the properties of a skewon-axion medium. 

The following alternative representation of the constitutive tensor is useful 
in many derivations and for a comparison with literature, see Post [50] , 

X m ■= \ ^ Kmn M , (32) 

with 

x ijkl = (1) X W + £ ym[fc pj] _ e klm[i p m j\ + ^3kl^ _ (33) 
20, principal 15, skewon 1, axion 

It is convenient to consider the excitation H and the field strength F as 6- 
vectors, each comprising a pair of two 3-vectors. The spacetime relation then 
reads 

7~ta \ _ f £ a 23fea \ ( ~Eb 
V 

Accordingly, the constitutive tensors are represented by the 6x6 matrices 



V a J I yba jj a J I ftb 



(34) 



K K _ f %>ba\ IK _ f^ab ® a* \ 

K I - orba mall X ~ \ (fa , qlab ■ l J0 J 



The 3x3 matrices 21, 58, £, D are defined by 

yba ;= ^OaOb } ;= 1 ^ h ^ x cdef _ ( g 6) 

£\ := \ hed X cMa , S)a b := \ eacd X° bcd , (37) 



or explicitly, recalling the irreducible decomposition (33) , 
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2l ab = -e ab - e abc ® a6 = ^ + e abc &> c , (38) 

€\ = 7 °6 - - 5 a b fi c c ) + adl (39) 

£a b = l\ +(Pa b -S b a fi c c ) + a5 b a . (40) 



The constituents of the principal part are the permittivity tensor e ab = e ba , 
the impermeability tensor \i~ h x = /x^ 1 , and the magnetoelectric cross-term 7°(,, 
with 7 C C = (Frcsnel-Fizeau effect). The skewon $ b a and the axion a describe 
electric and magnetic Faraday effects and (in the last two relations) optical 
activity. If we substitute (38), (39), (40) into (34), we find a 3-dimcnsional 
explicit form of our Axiom 5 formulated in (24): 

V a = (e ab - e abc fi c °) E b + ( 7 a b + fi b a - 5% fi c c ) B b + aB a , (41) 
= - e abc ^ C ) B b + (- 7 b a + $ a b - S b E b -aE a . (42) 

10 Propagation of electromagnetic rays ("light") 

After the spacetime relation (Axiom 5) has been formulated, we have a com- 
plete set of equations describing the electromagnetic field. We can now study 
the propagation of electromagnetic waves a la Hadamard. The sourceless 
Maxwell equations read 

dH = , dF = . (43) 

In the geometric optics approximation (cquivalcntly, in the Hadamard ap- 
proach) an electromagnetic wave is described by the propagation of a discon- 
tinuity of the electromagnetic field. The surface of discontinuity S is defined 
locally by a function & such that # = const on S. The jumps [ ] of the elec- 
tromagnetic quantities across S and the wave covector q := d<P then satisfy 
the geometric Hadamard conditions: 

[H]=0, [dH] = q A h = => h = q A c , (44) 
[F] =0, [dF] =qAf = f = qAa. (45) 

Here c and a are arbitrary 1-forms. 

We use the spacetime relation and find for the jumps of the field derivatives 

h = «(/) = «(/) + af , (46) 

with k := W/c + ( 2 )k. Accordingly, 5 

qAh = qAK,(qAa) = 0. (47) 

This equation is a 3-form with 4 components. We have to solve it with respect 
to a. As a first step, we have to remove the gauge freedom a — > a + qip present 

5 Compare the corresponding tensor analytical formula Opx" 13 "' 6 d-yAs = (see 
Post [50], Eq.(9.40) for x ^-<^ = 0). 



12 Friedrich W. Hehl and Yuri N. Obukhov 

in (47). We choose the gauge $° = q. After some heavy algebra, we find (see 
[20] for details, a, b, ... = 1,2,3) 

W ab a b = 0, with W ab := x 6a " b ■ (48) 
These are 3 equations for three a^sl Nontrivial solutions exist provided 

W := det W ab = ^ia bc e d efW ad W be W cf = 0. (49) 

We can rewrite the latter equation in a manifestly 4-dimensional covariant 
form (e abc = e 6abc , e;° = qi ), 

W = J e mnpg e rstu x m X jpsk X lqtu HWM = , 

with 9 := det(ej a ). The 4-dimensional tensorial transformation behavior is 
obvious. 

Wc define 4th-order Tamm-Rubilar (TR) tensor density of weight +1, 



m kl (-v\ — —e f , v mnr(i j\ps\k l)qtu 

v \XJ ■ — c mnpq c rsiu A A A 



(50) 



It is totally symmetric G l -* kl {x) = Q^ kl \x)- Thus, it has 35 independent 
components. Because x^ kl = X 1 ^ 1 + ae l i kl , the total antisymmetry of e yields 
Q(x) = Q{x)- An explicit calculation shows that 

G ijkl {x) = Q ijkl {^x) + (1) X m(llnU Pn ] ■ (51) 

Summarizing, we find that the wave propagation is governed by the ex- 
tended Fresnel equation that is generally covariant in 4 dimensions: 



g ijkl (x)qiqjq k qi=0. 



(52) 



The wave covectors q lie on a quartic Fresnel wave surface, not exactly what we 
are observing in vacuum at the present epoch of our universe. Some properties 
of the TR-tensor, see [53], were discussed recently by Beig [3]. 

Extended Fresnel equation decomposed into time and space 

Recalling the '6-vector' form of the spacetime relation (34) with the 3x3 
constitutive matrices (36) and (37), we can decompose the TR-tensor into 
time and space pieces: G omo =: M , g m0a =: \M a , g 0Qab =: \M ab , g 0abc =: 



\M abc , g abcd =: M abcd . Then the Fresnel equation (52) reads 



q* +ql q a M a +q 2 q a q b M ab +q q a q b q c M abc + q a q b q c q d M abcd = , (53) 
M o Mi M 2 M 3 M 4 
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or 



with 



M qt + Ah ql + M 2 q 2 a + M 3 q + M 4 = , 



'e I ' 



M = det 21 , M a = -hd (2l fca 2l ce € d e + 2l oh 2l ec £ 
+(£ d c + D c d )(^ a € b) d + T) d {a % b)c ) - € d d ^i c{a € b) c 

M abc = e de(c\ \<8 df (%ab) ®J _ Jjaypf ) + <B /d (2l ab) <£ f e - 21/ 1 Q £ b ;> ] 

£ e b) Dj + S/ 



(54) 
(55) 



y^abcd _ e ef(c £ \gh\d sg 



-B a6 > ^e-C^' 



(56) 

(57) 
(58) 




Fig. 4. Fresnel wave surface for anisotropic permittivity e ab = diag(39.7, 15.4, 2.3) 
and trivial impermeability fjT^ — /i ( ^ 1 diag(l, 1, 1). The skewon field vanishes. There 
are two branches, the outer part of the surface is cut into half in order to show 
the inner branch. We use the dimensionless variables x := cqi/qo, y := cqz/qo, z := 
cqz/qo- 
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Presnel wave surfaces 

Let us look at some Fresnel wave surfaces in order to get some feeling for 
the physics involved. Divide (53) by q$ (here q is the frequency of the wave) 
and introduce the dimensionless variables (c = velocity of light in special 
relativity) 

x a := c— . (59) 

Then we have 

M + x a h x a x b — s- + x a x b x c — 3 h x a x b x c x d — t — = . (60) 

c c 1 C A c 4 

We can draw these quartic surfaces in the dimensionless variables x = x\, 
y = X2, z = x 3 , provided the M's are given. According to (55)-(58), the M's 
can be expressed in terms of the 3x3 matrices 21, 05, £, !D. These matrices are 
specified in (38)-(40) in terms of the permittivity etc.. A comparison with the 
spacetime relations in the form of (41), (42) is particularly instructive. 



1 




Fig. 5. Fresnel wave surface for trivial permittivity e ab = eodiag(l, 1, 1) and trivial 
impermeability = ^i^ 1 diag(l, 1, 1) with a skewon field of electric Faraday type 
$ 3 ° — 3.1Ao (all other components vanish). The surface has the form of a toroid 
(depicted with two cuts). We use the dimensionless variables x := cqi/qo, y := 
cq 2 /qo, z := cq 3 /q a . 

Let us start with a simple example. We assume that the permittivity is 
anisotropic but still diagonal, e ab = diag(ei, £2, £3), whereas the imperme- 
ability is trivial fj,~ b = /i Q " 1 diag(l, 1, 1). No skewon field is assumed to exist. 
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Whether an axion field is present or not doesn't matter since the axion does 
not influence the light propagation in the geometrical optics limit. With Math- 
ematica programs written by Tertychniy [59] , we can construct for any values 
of ei, £2, £3 the Fresnel wave surface; an example is displayed in Fig. 4. 

More complicated cases are trivial permittivity e ab — £0 diag(l, 1, 1) and 
trivial impermeability = fj,Q diag(l, 1, 1), but a nontrivial skewon field. 
We can take a skewon field of electric Faraday type j$ 3 °, for example, see 
Fig. 5, or of magnetoelectric optical activity type jS-^ 2 = ^ 2 1 ^ see Fig. 6. In 
both figures and in the subsequent one Ao = -\/£o/a*o is the admittance of 
free space. The characteristic feature of the skewon field is the emergence of 
specific holes in the Fresnel surfaces that correspond to the directions in space 
along which the wave propagation is damped out completely [45] . This effect 
is in agreement with our earlier conclusion on the dissipative nature of the 
skewon field. 



Fig. 6. Fresnel wave surface for trivial permittivity e ab = eodiag(l, 1, 1) and trivial 
impermeability fj,~} — /i ( ^ 1 diag(l, 1, 1) with a skewon field of the magneto-electric 
optical activity type j5\ 2 =^2 1 = 0.8 Ao (all other components vanish). It has two 
intersecting toroidal branches. We use the dimensionless variables x := cqi/qo, y := 
cq 2 /qo, z := cq :i /q . 

Now we can combine anisotropic permittivity with the presence of a skewon 
field. Then we expect to find some kind of Fig. 4 "enriched" with holes induced 
by the skewon field. This time we choose a spatially isotropic skewon field with 
#1 =y g- 2 2 =$5 3 = _l ^ q o ^ 0. The outcome is depicted in Fig.7. The four 
holes confirm our expectation. 






Fig. 7. Fresnel wave surface for anisotropic permittivity s ab = diag(2.4, 14.8, 54) 
and trivial impermeability fi~£ = |zjj" 1 diag(l, 1, 1) with a spatially isotropic skewon 
field fii 1 =fi 2 2 =^3 3 = -| = 0.25 A (all other components vanish). We use 
the dimensionless variables x ~ cqi/qo, y := cq^jqo,, z := cq^/qo. 



11 No birefringence in vacuum and the light cone 

The propagation of light in local and linear premetric vacuum electrodynamics 
is characterized by the extended Fresnel equation (52) or (54). We can solve 
the Fresnel equation with respect to the frequency q , keeping the 3-covector 
q a fixed. With the help of Mathematica, we found the following four solutions 
[33]: 



= y/a±Jp + -l=-8, (61) 



We introduced the abbreviations 
if a 

a := 



12M ( 



T + (b + y/c) 3 - 2M 2 )+S 2 , (63) 



V (b + V$ 

P--=T^T { a — - ( b + Vc) i - 4M 2 j + 2S 2 , (64) 

' (2SM 2 - M 3 ) - 2S 3 , S := , (65) 



4M v ' 4M 
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with 



a := 12M M 4 - 3MiM 3 
27 
2 

c := 4 (6 2 - a 3 ) . 



(66) 



b := — M M| - 36M M 2 M 4 - -M 1 M 2 M 3 + —M\M^ + M 2 3 , (67) 

(68) 



Vanishing birefringence 

Now, let us demand the absence of birefringence (also called double refrac- 
tion). 6 In technical terms this means, see the solutions (61), (62), that j3 = 
and 7 = 0. Then we have the degenerate solution 

t i— Mi l Mr 

%= ^-4Mo' 9 ° = ^"4Mo- (69) 

The condition 7 = yields directly M 3 = Mi (4M M 2 - M 2 ) /8M 2 , and, 
using this, we find 

3M 2 - 8M M 2 



Thus, 



iTi _ + / 3M 2 - 8M M 2 Mi 
16M 2 4Mo 

Accordingly, the quartic wave surface (54) in this case reduces to 

[(«b-«5)(«b-9o)] 2 = 0- (72) 

Multiplication yields 

2 , 1 Mi , 1 M 2 1 f Mi N J 



^2^* + ^-8UJ = °' (?3) 

If we substitute M ,Mi,M 2 as defined in (53), we have explicitly = 
0,1,2,3) 

, 1 M a If M ab M a M b \ , . 

9 v m ■= <& + 2 w »«. + s\ ~W ~ ) qa% = ' (74) 

This equation is quadratic in the 4-dimcnsional wave covector Therefore 
we recover the conventional light cone of general relativity at each point of 
spacetime, see Fig. 8. Thereby the causal structure of spacetime is determined. 



Similar considerations on vanishing birefringence, for weak gravitational fields, 
are due to Ni [44]. He was also the first to understand that the axion field doesn't 
influence light propagation in the geometrical optics limit. 
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Thus, up to a scalar factor, we derived the Riemannian metric of general 
relativity. 

Moreover, as we have shown [20, 25], we find the correct Lorentzian sig- 
nature. The Lorentzian (also known as Minkowskian) signature can be traced 
back to the Lenz rule, which determines the sign of the B term in the induc- 
tion law. 7 And this sign is different from the one in the corresponding T> term 
in the Oersted-Ampere-Maxwell law. In other words, the Lorentz signature 
is encoded in the decomposition formulas (9) and (14). Neither is the minus 
sign in (9) a convention nor the plus sign of the E A da term in (14). Since 
the Lenz rule is related to the positivity of the electromagnetic energy, the 
same is true for the Lorentzian signature. This derivation of the signature of 
the metric of spacetime from electrodynamics provides new insight into the 
structures underlying special as well as general relativity. 



Fig. 8. Null cones fitted together to form a conformal manifold (see Pirani and 
Schild [49]). 



Usually it is argued that the signature should be derived from quantum field 
theoretical principles; for a corresponding model, see, e.g., Froggatt & Nielsen [15]. 
Needless to say that it is our view that classical premetric electrodynamics together 
with the Lenz rule and a local and linear spacetime relation is all what is really 
needed. 
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12 Dilaton, metric, axion 

At first the skewon and the axion emerged at the premetric level in our theory 
and only subsequently the metric. Consequently, the axion and the skewon 
should be regarded as more fundamental fields (if they exist) than the metric. 
In the meantime, we phased out the skewon field since we insisted, in Sec.ll, 
on vanishing birefringence in vacuum. 

As to the metric, we recognize that multiplication of the metric by an 
arbitrary function X(x) was left open in the derivation of the last section, see 
(74): 

\(x)g i *(x)q i q j =0. (75) 

Thus, only the conformally invariant part of the metric is determined. In 
other words, we have actually constructed the conformal (or the light cone) 
structure on the spacetime manifold, see, e.g., Weyl [63, 64], Schouten [55], 
and Pirani & Schild [49]. 

It is known from special relativity that the light cone (with Lorentzian 
signature) is invariant under the 15-parameter conformal group, see Barut & 
Raczka [2] and Blagojevic [4]. The latter, in Minkowskian coordinates x l , is 
generated by the following four sets of spacetime transformations: 

Translations (4 param.) x l — > x % = x % + a 1 , (76) 
Lorentz transf. (6 param.) x l — > x l = A 1 j , (77) 
dila(ta)tion (1 param.) x % — > x 1 = px\ (78) 

prop. conf. transf. (4 param.) x l — > x l = : : — -. (79) 

1 -|- X 3 ~\~ KjK^ X 

Here a l ,A l j 7 p, n l are the 15 constant parameters, and x 2 := gijX l x 3 . The 
Poincare subgroup (76), (77) (for a modern presentation of it, see Giulini 
[17]) leaves the spacetime interval ds 2 = g i jdx l dx : > invariant, whereas the 
dilatations (78) and the proper conformal transformations (79) change the 
spacetime interval by a scaling factor ds 2 — > p 2 ds 2 and ds 2 — > <j 2 ds 2 , respec- 
tively (with cr _1 := 1 + 2kj .t j + k^k- 7 x 2 ). In all cases the light cone ds 2 = 
is left invariant. The Weyl subgroup, which is generated by the transforma- 
tions (76)- (78), and its corresponding Noether currents were discussed by, e.g., 
Kopczyhski et al. [30]. 

For massless particles, instead of the Poincare group, the conformal or 
the Weyl group come under consideration, since massless particles move on 
the light cone. Even though the light cone stays invariant under all trans- 
formations (76)-(79), two reference frames that are linked to each other by a 
proper conformal transformation don't stay inertial frames since their relative 
velocity is not constant. If one wants to uphold the inertial character of the 
reference frames, one has to turn to the Weyl transformation, that is, one has 
to specialize to n l = 0. 

The conformal group in Minkowski space illustrates the importance of the 
light cone structure on a flat manifold. This is suggestive for the light cone on 
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an arbitrarily curved manifold, even though there is no direct relation between 
(76)-(79) and the light cone structure we derived in the last section. 

The light cone metric g^ introduces the Hodge star * operator. We then 
can straightforwardly verify that the principal part of the spacetime relation 
is determined as H ~ *F, where the coefficient of proportionality can be an 
arbitrary scalar function A(x) of the spacetime coordinates. This function is 
naturally identified with the dilaton field, see Brans [6] and Fujii & Maeda 
[16]. Introducing the (Levi-Civita) dual of the excitation, H % i := | e n kl H^i, 
we can then finally rewrite the spacetime relation for vanishing birefringence 
in vacuum as 

fri - [ V=gg tk (x) g 3l (x) + ^ z m ] F kl , (80) 

dilaton axion 

that is, we are left with the constitutive fields dilaton A, metric g lJ , and axion 
a. The combination yj— g g^ k (x) g l ^(x) is conformally invariant, in complete 
agreement with the above analysis. 



13 Setting the scale 

The conformal structure of spacetime is laid down in (74) . Hence only 9 of the 
10 independent components of the pseudo-Riemannian metric gij are specified. 
We need, in addition to the conformal structure, a volume measure for arriving 
at a unique Riemannian metric. This can be achieved by postulating a time 
or length standard. 

In exterior calculus, (80) reads 

H = X(x) *F + a(x) F . (81) 

The axion has not been found so far, so we can put a — 0. Moreover, under 
normal cicumstances, the dilaton seems to be a constant field and thereby 
sets a certain scale, i.e., X(x) = Ao, where Ao is the admittance of free space 8 
the value of which is, in Si-units, 1/(377 fi). (The exact implementation of 
this assumption will have to be worked out in future.) Accordingly, we are left 
with the spacetime relation of conventional Maxwell-Lorentz electrodynamics 



H = X *F or it* = \ sf=g~g ik {x)g> l {x)F kl =A„^i™. (82) 



14 Discussion 

Weyl [63, 64], in 1921, proved a theorem that the projective and the conformal 
structures of a metrical space determine its metric uniquely. As a consequence 

8 Our electrodynamical formalism is independent of the chosen system of units, 
as we discussed elsewhere [22]. 
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Weyl [63] argued that ...the metric of the world can be determined merely by 
observing the "natural" motion of material particles and the propagation of 
action, in particular that of light; measuring rods and clocks are not required 
for that. Here we find the two elementary notions for the determination of the 
metric: The paths of a freely falling point particles, yielding the projective 
structure, and light rays, yielding the conformal structure of spacetime. Later, 
in 1966, Pirani and Schild [49], amongst others, deepened the insight into the 
conformal structure and the Weyl tensor. 

In 1972, on the basis of Weyl's two primitive elements, Ehlers, Pirani, and 
Schild (EPS) [12] proposed an axiomatic framework in which Weyl's concepts 
of free particles and of light rays were taken as elementary notions that are 
linked to each other by plausible axioms. Requiring compatibility between 
the emerging projective and conformal structures, they ended up with a Weyl 
spacetime 9 (Riemannian metric with an additional Weyl covector) . They set a 
scale [as we did in the last section, too] and arrived at the pseudo- Riemannian 
metric of general relativity. In this sense, EPS were able to reconstruct the 
metric of general relativity. 

Subsequently, many authors improved and discussed the EPS-axiomatics. 
Access to the corresponding literature can be found via the book of Majer 
and Schmidt [40] or the work of Perlick [47, 48] and Lammerzahl [31], e.g.. 
For a general review one should compare Schelb [54] and for a new axiomatic 
scheme Schroter [56]. 

As stated, the point particles and light rays were primary elements that 
were assumed to exist and no link to mechanics nor to electrodynamics was 
specified. The particle concept within the EPS-axiomatics lost credibility 
when during the emergence of gauge theories of gravity (which started in 
1956 with Utiyama [61] even before the EPS-framework had been set up in 
1972) the first quantized wave function \P for matter entered the scene as 
an elementary and "irreducible" concept in gravity theory. When neutron in- 
terference in an external gravitational field was discovered experimentally in 
1975 by Collella, Overhauser, and Werner (COW) [7], see also [52], Sec.7, it 
was clear that the point particle concept in the EPS-framework became un- 
tenable from a physical point of view. For completeness let us mention some 
more recent experiments on matter waves in the gravitational field or in a 
noninertial frame: 

• The Werner, Staudcnmann, and Colella experiment [62] in 1979 on the 
phase shift of neutron waves induced by the rotation of the Earth (Sagnac- 
type effect), 

• the Bonse & Wroblewski experiment [5] in 1984 on neutron intcrferometry 
in a noninertial frame (verifying, together with the COW experiment, the 
equivalence principle for neutron waves), 



9 Time measurement in Weyl spacetime were discussed by Perlick [47] and by 
Teyssandier & Tucker [60]. 
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• the Kasevich & Chu interferometric experiment [27] in 1991 with laser- 
cooled wave packets of sodium atoms in the gravitational field, 

• the Mewes et al. experiment [41] in 1997 with interfering freely falling 
Bose-Einstein condensed sodium atoms, see Ketterle [28], Fig. 14 and the 
corresponding text, 

• the Nesvizhevsky et al. experiment [42, 43] in 2002 on the quantum states 
of neutrons in the Earth's gravitational field, and 

• the Fray, Hansch, et al. [14] experiment in 2004 with a matter wave in- 
terferometer based on the diffraction of atoms from effective absorption 
gratings of light. This interferometer was used for two stable isotopes of 
the rubidium atom in the gravitational field of the Earth. Thereby the 
equivalence principle was tested successfully on the atomic level. 

Clearly, without the Schrodinger equation in an external gravitational field 
or in a nonincrtial frame all these experiments cannot be described. 10 Still, in 
most textbooks on gravity, these experiments are not even mentioned! 

In the 1980's, as a reaction to the COW-experiment, Audretsch and 
Lammerzahl, for a review see [1], started to develop an axiomatic scheme for 
spacetimc in which the point particle was substituted by a matter wave func- 
tion and the light ray be a wave equation for electromagnetic disturbances. 
In this way, they could also include projective structures with an asymmetric 
connection (i.e., with torsion), which was excluded in the EPS approach a 
priori. 

Turning to the conformal structure, which is in the center of our interest 
here, Lammerzahl et al. [32], see also [51, 18], reconsidered the Audrctsch- 
Lammcrzahl scheme and derived the inhomogeneous Maxwell equation from 
the following requirements: a well-posed Cauchy problem, the superposition 
principle, a finite propagation speed, and the absence of birefringence in vac- 
uum. The homogeneous Maxwell equation they got by a suitable definition 
of the electromagnetic field strength. With a geometric optics approximation, 
compare our Sec. 10, they recover the light ray in lowest order. And this is the 
message of this type of axiomatics: Within the axiomatic system of Audretsch 
and Lammerzahl et al., the light ray, which is elementary in the EPS-approach, 
can be derived from reasonable axioms about the propagation of electromag- 
netic disturbances. As with the substitution of the mass point by a matter 
wave, this inquiry into the physical nature of the light ray and the correspond- 
ing reshaping of the EPS-scheme seems to lead to a better understanding of 
the metric of spacetimc. And this is exactly where our framework fits in: We 
also build up the Maxwell equations in an axiomatic way and are even led to 
the signature of the metric, an achievement that needs still to be evaluated in 
all details. 



A systematic procedure of deriving the COW result by applying the equivalence 
principle to the Dirac equation can be found in [19]. 
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15 Summary 

Let us then summarize our findings: We outlined our axiomatic approach to 
electrodynamics and to the derivation of the light cone. In particular, with 
the help of a local and linear spacetime relation, 

• we found the skewon field (15 components) and the axion field a 
(1 component), 

• we found a quartic Fresnel wave surface for light propagation. 

• In the case of vanishing birefringence, the Fresnel wave surface degener- 
ates and we recovered the light cone (determining 9 components of the 
metric tensor) and, together with it, the conformal and causal structure 
of spacetime and the Hodge star * operator. 

• If additionally the dilaton A (1 component) is put to a constant, namely 
to the admittance of free space A [1/(377 fi) in SI- units], and the axion a 
removed, we recover the conventional Maxwell-Lorentz spacetime relation 
H = X *F. 

Thus, in our framework, the conformal part of the metric emerges from the 
local and linear spacetime relation as an electromagnetic construct. In this 
sense, the light cone is a derived concept. 
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